library(tidyverse)

## Read Nexis data
  
out_full <- read_rds("data/nexis.rds")

## Agg for plot

plot_df <- out_full %>% 
  group_by(year) %>% 
  summarise(count = n())

# Figure A.26 - Newspaper articles ----

p1 <- ggplot(plot_df, aes(year, count)) +
  geom_bar(stat = 'identity', width = 0.65, 
           fill = 'grey95', color = 'black') +
  scale_x_continuous(breaks = unique(plot_df$year)) +
  ylab("No. of newspaper articles covering\nthe new census methods") +
  xlab('') +
  theme_bw() +
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1))
p1